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Abstract. We prove a uniqueness result of solutions for a system of 
PDEs of Monge-Kantorovich type arising in problems of mass transfer 
theory. The results are obtained under very mild regularity assumptions 
both on the reference set C R™, and on the (possibly asymmetric) 
norm defined in Q. In the special case when H is endowed with the 
Euclidean metric, our results provide a complete description of the sta- 
tionary solutions to the tray table problem in granular matter theory. 



1. Introduction 

The model system usually considered for the description of the stationary 
configurations of sandpiles on a tray table is the Monge-Kantorovich type 
system of PDEs 

— div(f Du) = f in SI, 
\Du\ < 1, v > in O, 
(l-\Du\)v = Q in ft, 
^u = <j> on dVt 

(see, e.g., @, |lO|, ^3|). The data of the problem are the flat surface of the 
table J!CR 2 , the profile of the tray table and the density of the source 
/ > 0. The dynamical behavior of the granular matter is pictured by the 
pair (u, v), where u is the profile of the standing layer, whose slope has not 
to exceed a critical value (\Du\ < 1) in order to prevent avalanches, while 
v > is the thickness of the rolling layer (see Figure ||). The condition 
(1 — \Du\)v = corresponds to require that the matter runs down only in 
the region where the slope of the heaps is maximal. On the border of the 
table the standing layer fills the gap with the height <f>, and the exceeding 
sand falls down. If 4> = 0, the model reduces to the open table problem, 
largely investigated in recent years (see e.g. 0, |U|, H, and the references 
therein) . 
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Figure 1 . The tray table problem 



In the present paper we consider a more general system of PDEs in the 
open, bounded and connected set O C W 1 



(1) 



- div(u Dp(Du)) = f inO, 

p(Du) <l,v>0 in 0, 

(l-p(Du))v = inO, 

u = 6 on dO, . 



in the unknowns v £ L\(0), u £ W 1 ' 00 ^) n C(0). The data are: 

- K C W l , a compact convex set of class C 1 containing the origin in its 
interior, and its (strictly convex) polar set K°. The set K appears in 
(Q) by means of its gauge function p: W 1 — > [0, +oo). Moreover the 
closure of O will be equipped with a (possibly asymmetric) geodesic 
distance cIl given in terms of the gauge function p° of the polar set 
K° (see Section g); 

- / G L\(n); 

- (ft: 80, — > M n , 1-Lipschitz function w.r.t. <1l. 

This general version takes into account the possibility of homogeneous aniso- 
tropics, and can be applied, for example, to the Bean's model for the de- 
scription of the macroscopic electrodynamics of hard superconductors (see 
[16, 17, |l^, [DJ, for the case cj> = and O, K smooth sets), and to optimal 
shape design (see M). 

S. Bianchini in |H] proved that the Lax-Hopf function (defined in (O) 
below), coupled with an explicit vj G L\(0) solves (|l|) (see Theorem |6,9| ). 
The function Vf is constructed as follows. The set O can be covered by a 
family of disjoint transport rays, and every f-component of a solution to (|l|) 
satisfies a first order linear ODE along almost every ray. The function vj 
is then uniquely defined along each ray, by requiring that it vanishes at the 
final point. This procedure gives a solution to (|l|), thanks to a disintegration 
formula for the Lebesgue measure along the rays. 

In this paper we shall provide a complete characterization of the solutions 
(u, v ) to ([[]). Starting from the fact that (u^,,Vf) is a solution, we show that 
a pair (n, v) is a solution if and only if both (u,vj) and {u^, v) are solutions, 
so that we can decouple the original problem (see Section ||). Moreover, the 
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admissible u-components are characterized as the solutions of a minimiza- 
tion problem which has an explicit minimal solution uj depending on the 
source term /. Since is the maximal solution of the problem, the unique- 
ness of the u-component corresponds to the case when Uf = u^. We prove 
that this can occur if and only if the final points in Q of the transport rays 
are contained in the support of the source / (see Section [?]). Concerning 
the f-component, in || it is proved that every admissible v has to vanish at 
the final points in f2 of the transport rays. Hence, by the very definition of 
Vf, it is the unique u-component if the union T of the rays with both end- 
points on <9f2 has null Lebesgue measure. We shall prove that the converse 
is also true: if the set T has positive Lebesgue measure, we can construct a 
nontrivial function v E L+(Q) such that — div(v Dp{Du r f ) )) = 0, so that the 
uniqueness of the v-component fails. 

These results generalize those in [16, FT?], Il8| to the case when the ref- 



erence set f2 is possibly non regular, the boundary datum cj) is possibly 
non- homogeneous, and the constraint K is a possibly non-strictly convex 
set. In this case we have to use a completely different strategy with respect 
to the previous approaches. 

Namely, if Q and K are smooth enough, then it is possible to define a 
suitable notion of (anisotropic) curvatures of d£l, and the analysis of solu- 
tions to d^can be performed using the methods of Differential Geometry 
(see, e.g., [13, |l6[ [H], |J, [25|, |27|]). In the non-smooth case these methods 



do not work, and a different approach is needed, based on the techniques 
developed in Optimal Transport Theory (see ||, |3l]]). Namely, system (|l]) 
can be understood as the PDE formulation of an optimal transport problem 
with strictly convex cost function p° , assigned initial distribution \i = / dx 
and transport potential given by the Lax-Hopf function ua, (so that the fi- 
nal mass distribution v is the measure concentrated on d£l induced by the 
transportation of p along the rays associated to u^). From this point of 
view, the minimization problem considered in Section|5] turns out to be the 
dual formulation of the Monge problem (see, e.g., |HJ for an exhaustive 
presentation of this subject). 

A comment on the regularity assumption on K is in order. Namely, if 
K is of class C , then its polar set K° need not be of class C , but it 
is strictly convex. It is well-known that the strict convexity of K°, which 
is the unit ball of the dual norm p°, allows non-trivial simplifications in 
the proofs; for example, the direction of the transport rays coincides with 
Dp(Du) almost everywhere in £1. Related results in the case of a general 
possibly non-regular constraint K can be found in Q. 

The paper is organized as follows. In Section [2] we collect all the nota- 
tion and preliminary results. In Section || we recall some generalization of 
the McShane's extension theorem for Lipschitz functions with constrained 
gradient. These general results will be used for introducing the maximal 
and the minimal extension of the boundary datum <j> consistent with the 
gradient constraint. Section [|] is devoted to the statement of the problem, 
some remarks on the hypotheses, and the statement of the existence result 
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due to S. Bianchini (see |], Thm. 7.1]). One of the main tools needed for 
proving the uniqueness results is the fact that ((!]) is the Euler-Lagrange 
necessary condition for a minimization problem with a gradient constraint. 
This property is proved in Section ||. The necessary and sufficient condi- 
tions for the uniqueness of the solutions (u, v) to (|l]) are proved in Sections ^ 
and ^ As a byproduct of our analysis, with a little additional effort, we 
provide a complete description of the singular sets considered in this kind 
of optimal transport problems. 

Acknowledgements. The authors wish to thank the anonymous referee 
for valuable comments that have improved the presentation of the manu- 
script. 

2. Notation and preliminaries 

The standard scalar product of x, y 6 M. n will be denoted by (x, y), while 
\x\ will denote the Euclidean norm of x. Concerning the segment joining x 
with y, we set 

fx, y] := {tx + (1 - t)y; t G [0, 1]}, jx, y{ := fx, yj \ {x, y}. 

Given a set A C M n , its interior, its closure and its boundary will be 
denoted by int A, A and dA respectively. 

A bounded open set O C W l (or, equivalently, O or dO) is of class C k , 
k G N, if for every point x$ G dO there exists a ball B = B r (xo) and 
a one-to-one mapping ij): B — > D such that ip G C k (B), if)^ 1 G C k (D), 

ip{B n O) c {x g W 1 - x n > o}, ip(B n dO) c{ie M n ; Xn = 0}. 

Let us now fix the notation and the basic results concerning the convex 
set which plays the role of gradient constraint for the u-component in (jl|). 
In the following we shall assume that 

(2) K is a compact, convex subset of W 1 of class C 1 , with E int K. 

Let us denote by K° the polar set of K, that is 

K° := {p G R n ; (p, x) < 1 Vx G K} . 

We recall that, if if satisfies @, then ^° is a compact, strictly convex 
subset of W 1 containing the origin in its interior, and K 00 = (K°)° = K 
(see, e.g., @). 

The gauge function p : R n — > R of if is defined by 

p{0 := M{t > 0; £etK} = max{(e, r?), r? G K } , (eR n . 

It is straightforward to see that p is a positively 1-homogeneous convex 
function such that K = {£ G M n : p(£) < 1}. The gauge function of the set 
K° will be denoted by p°. 

The properties of the gauge functions needed in the paper are collected 
in the following theorem. 
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Theorem 2.1. Assume that K C W 1 satisfies (^). Then the following hold: 

(i) p is continuously differ entiable in R n \ {0} ; and 

p(Z + v)<p(0 + p(v) v^er. 

(ii) K° is strictly convex, and 

P°^ + V) = P°iO+ P°(V) & 3A>0: £ = Xn or n = 
(Hi) For every £ 7^ 0, Dp{^) belongs to dK° , and 

(Dp(0, = p(0, (p, < p(0 Vp eK°, P ^ D P {i) . 

Proof. See M, Section 1.7. □ 



In what follows we shall consider M. n endowed with the possibly asym- 
metric norm p°(x — y), x,y G W 1 . By Theorem 2.1(h), the unit ball K° 



of p° is strictly convex but, under the sole assumption (0), it need not be 
differentiable. Moreover, the Minkowski structure (M. n ,p ) is not a metric 
space in the usual sense, since p° need not be symmetric (for an introduc- 
tion to nonsymmetric metrics see [22|). Finally, since K° is compact and 



G mtK°, then the convex metric is equivalent to the Euclidean one, that 
is there exist c\, C2 > such that ci|£| < p°(C) < c 2|£| for every £ G M n . 

The main motivation for introducing the convex metric associated to p° is 
the fact that the Sobolev functions with the gradient constrained to belong 
to K are the locally 1-Lipschitz functions with respect to p°, as stated in 
the following result. 

Lemma 2.2. Let O C W 1 be a nonempty open set, and assume that the 
set K C W 1 satisfies (^). Let p° be the gauge function of K° . Then the 
following properties are equivalent. 

i) u: O — > R is a locally 1-Lipschitz function with respect to p°, i.e. 
(3) ^(2:2) — u(x\) < p°(x2 — x\) for every [xi,^] C O. 

ii) u G W^°°(0), and Du(x) G K for a.e. xeO. 

Proof. Let u be a locally 1-Lipschitz function in O with respect to the metric 
p°. Then u is locally Lipschitz in 0, and if xq is a differentiability point of 
u, by (|3|) we have 

o(|*|) = ^o + *0-"(^o)-*W^).0 < p o {0 _ {Du{xq)j 

< 1 - (Du(x ), 0, V£ G K°. 

Hence {Du(x ), f) < 1 for all £ G K°, that is Du(x ) G (K°)° = K. Since 
u is a locally Lipschitz function in O with gradient in the bounded set K, 
(ii) holds. 

Conversely, let u satisfy (ii), and let [xi,^] C O. Clearly, we can assume 
x\ 7^ X2- Let S r := {z G B r (0) : (z — x\, X2 — x\) = 0}, and let Q r := 
{U l x i + z,X2 + z\ : z G S r } be the cylinder of radius r around the segment 
[xi,X2]. Let us choose r > such that Q r C By (ii) there exists 
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a set N C S r with (n — 1)-Lebesgue measure zero, such that for every 
z G S r \ N, u is absolutely continuous along the segment {xi + z,x 2 + zj, 
and Du(x\ + z + s(x2 — x±)) G K for a.e. s G [0, 1]. Let (zj) C S r \ N, 
Zj — > 0. Then, for every j G N, 

it(a;2+^j)— u{x\+Zj) = / (D«(xi+^j+s(aj 2 — a^i)), »2— x{) ds < p°{x 2 -xi). 

Jo 

The conclusion now follows from the continuity of u in O, passing to the 
limit in j. □ 

In the sequel we need to consider locally 1-Lipschitz functions with re- 
spect to p° that agree with a given continuous boundary datum. This can be 
done only if the boundary values are compatible with the gradient bounds, 
and, in absence of assumptions on the regularity of the reference set, the 
compatibility condition concerns the variations of the data along every ad- 
missible path. 

Here we introduce some notation for the curves, while the compatibility 
condition will be discussed in Section ||. 

Definition 2.3. Let D C W 1 be a nonempty compact set. For every x\, x 2 G 
D, T Xl X2 will denote the (possibly empty) family of absolutely continuous 
curves 7: [0, 1] — >■ D such that 7(0) = x\ and 7(1) = x 2 - 

For every absolutely continuous curve 7: [0,1] —> M n , let us denote by 
L(7) its length with respect to the convex metric associated to p°, that is 

L( 7 ) := C P \l\t))dt. 
Jo 

If D is a compact subset of M. n , by a standard compactness argument we 
have that if X2 / 0, then there exists a geodesic in D joining x\ to x 2 , i.e. 
there exists a curve 7 G T Xl X2 such that £(7) < £(7) for every 7 G T X1 X2 



(see e.g. |, Thm. 4.3.2], |5|, §14.1]). 

For sake of completeness we briefly recall that any geodesic curve lying 
in the interior of D is in fact a segment, due to strict convexity of p°. 

Lemma 2.4. Let D C M™ be a nonempty compact set. Let xi,x 2 G D be two 
points such that F X1 X2 7^ ; and let 7 G F® l X2 be a geodesic. If 7(i) G int D 
for every t G (0, 1), then the support of ^ is the segment \x\,x 2 \. 

Proof. By an approximation argument we can assume that x\,x 2 G vcAD 
so that 7(t) G int-D for every t G [0, 1]. Let n(t) = x\ + t(x 2 — x\), t G [0, 1], 
be a parametrization of the segment \x\, x 2 J. Assume by contradiction that 
the support of 7 is not the segment jxi, x 2 J. Since K° is strictly convex we 
have that £(7) > L(jf) = p°(x 2 — xi). Since the support of 7 is contained 
in int-D, there exists e G (0, 1) such that the support of 7 e := 7 + e(r] — 7) 
is contained in intD. Then 

L{le)= / 1 p°((l-e)7 / + ^')< l\(l-e)p\ 1 ') + ep\r 1 ')] 
JO Jo 

= (l-e)L{ 1 )+eL(r } )<L( rf ), 
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in contradiction with the minimality of 7. □ 

3. LlPSCHITZ EXTENSIONS 

In this section, C and D will be two sets such that 

D is a nonempty compact subset of R n , and 
^ ^ C is a closed subset of D containing dD, 

so that D \ C is an open (possibly empty) set. 

For a given continuous function ip : C — > R, we want to discuss the exis- 
tence of a locally Lipschitz extension of ip in D with the gradient constrained 
to belong to a convex set K. 

It is well known that, in the case D = $7 bounded open subset of W 1 , 
C = (9f2, and K satisfying (|2|), a necessary and sufficient condition for the 
existence of such an extension is 

iP(x) - ip(y) < L( 7 ) Vx, y G V 7 € r", 

where Ty X is the family of absolutely continuous curves 7 : [0, 1] — > f2 such 
that j(t) G for every t G (0, 1), 7(0) = y and 7(1) = x (see [||, Sect. 5] 
and 0). Moreover, if has a Lipschitz boundary, one can relax the above 
necessary and sufficient condition by taking paths in (see p6|). 

In order to handle more general situations, we will consider a continuous 
function ip : C — > M satisfying 

(5) iP{x) - iP(y) < L( 7 ) Vx, y G C, V 7 G , 

where r£ ^ is the family of curves introduced in Definition 2.3 



This condition can be rephrased in the setting of the length space (D, di), 
where di '■ D x D — > [0, +00] is the (possibly asymmetric) distance defined 
by 

(6) d L (x, y) := inf{L( 7 ) : 7 G T° x }, x,y G D. 

(We refer to Chap. 2] for the basic properties of length spaces; see in 
particular Remark 2.2.6 for the extension to possibly asymmetric metrics.) 
Namely, (|5|) is equivalent to requiring that ip is a 1-Lipschitz function with 
respect to the d^. 

In the following lemma we recall the well known fact that (||) turns out 
to be a sufficient condition for the existence of minimal and maximal 1- 
Lipschitz extensions (w.r.t. di) of ip . Notice that (||) is no more a necessary 
condition, as it can be seen when D is a closed segment in M. 2 . 

Lemma 3.1. Let C, D satisfy (Oj, let K satisfy and let ip: C — > M. be 
a continuous function satisfying (^). Then the functions 

u*(x) := inf ty(y) + L( 7 ) : y G C, 7 G ^y,x} , x £ D, 
w*(x) := sup {iP{y) - L( 7 ) : y G C, 7 G 1^} , x £ D, 

are continuous in D, and = = ip in C . Moreover, u^,w^ G 
W^intZ?) and Du^^Dw^ £ K a.e. in intD. 
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As a consequence of Lemma 3J., the functions vft and belong to the 
functional space 

Z,p := {u G C(D) n W^QxA D): u = ip on C, Du G K a.e. in int D) . 

It should be noted that in the definition of u^(x) the infimum is taken over 
all paths in D joining y G C to x, while the supremum in the definition 
of is taken over all paths in D joining x to y G C . This asymmetry is 
needed in order to construct exactly the maximal and the minimal elements 
of Z^, as it is stated in Lemma below (see |], Prop. 2.1]). 

Definition 3.2. Given x G D, we shall say that y G C is a visible point 
for x if the segment Jy, x\ is entirely contained in the open set D\C. We 
shall denote by V£ the set of all visible points for x, that is V% = {y G C : 
}y,x{cD\C}. 

Lemma 3.3. The functions ift , vft defined in ^) are characterized in the 
following way: 

v?(x) = min{>(y) + p°{x-y) : y G V? }, x G D, 
w^[x) = max{^(y) - p°(y - x) : y G V^F}, x G D. 

Moreover vft < u < for every u G Zj, . 

Corollary 3.4. Let x G D \ C be a point of differentiability for (resp. 
w*). Then Du^(x) (resp. Dw^{x)) belongs to dK. 

Proof. Given x G D \ C, let y G V% be such that 

u^(x)=iP(y)+p°(x-y). 

Setting £ = (x — y)/p°(x — y), the fact that G Z$ and Lemma 2J1 imply 
that 

u^{x-si)-u^{y) <p°{x-y)-s, (x)-^ (x-s£) <s Vs G [0, p°(x-y)] . 
Then we have 

-s <u^(x- s£) - u^(x) = u^(x - sQ - u^{y) - p°(x -y)<-s, 
and hence 

(9) u^{x - sQ = u^{x) - s, Vs G [0, p°(x — y)] . 

If is differentiable at x, (||) implies that (Du^(x), £) = 1, so that 
p(Dv?(x)) = max{Du f (x), 77) > 1 . 

On the other hand, by Lemma 3.1, Du^(x) G K, hence Du^(x) belongs to 
dK. 

The proof for is omitted, since it is entirely similar. □ 
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4. The existence result 

Our main goal is to provide a full analysis of the Monge-Kantorovich 
system of PDEs 



(10) 



'- div(v Dp(Du)) = f inO, 

p(Du) < 1, v > a.e. in O, 

(1 — p(Du))v = a.e. in O, 

u = 4> pointwise on 90 , 

under the following assumptions: 

(HI) O is a nonempty bounded connected open subset of R n , n > 1. 
(H2) / G L 1 ($7), / > a.e. in O. 

(H3) p: R n -)• R is the gauge function of a set if C R n satisfying (||). 
(H4) </> : 90 — > R is a continuous function, satisfying the condition 

0(aO " < Mt) := f p°(j'(t)) dt Vx, y G 90, V 7 G T^, 

J o 

where 

Fy,* := = { 7 G AC7([0, 1],0) : 7 (0) = y, 7 (1) = x} . 

Remark 4.1. In the first equation in ( |l0"| ) we understand -Dp defined also at 
the origin since, by the third condition, v = a.e. where -Du = 0. 

The functional setting for the unknowns (u, v) in (0) is I^xi}(Q), where 

(11) := {u G C(O) n W 1,oc («) : u = on 90, Du £ K a.e. in 0} 
and 

L^(O) := {v G L x (0) : u > a.e. in 0} . 
By a solution of ( |To|) we mean a pair (u, u) G x Li_(0) satisfying 

(12) / v(Dp(Du), Dtp) = f ftp yip G C °°(O) 

and (1 — p(Du))v = a.e. in 0. 

Let us define the Lax-Hopf function ii,*:fi->Kby 

(13) u^x) := inf {(j){y) + L( 7 ) : y G 50, 7 G r„, x } , x G O, 

which corresponds to the function defined in ^ in the special case 
D = O, C = 90, and ip = <ft (observe that C = 90 D 90 = 3D and 
O C intD). 



As a consequence of Lemmas |3.1| , |3.3| , and Corollary |3.4| , we obtain the 
following result, which is well-known when O is a Lipschitz domain (see e.g. 
|2(| Chap. 5]). 

Theorem 4.2. If (HI), (H3) and (H4) are fulfilled, then G and 
Du^ G dK a.e. in 0. Moreover, is the maximal element in X^, and it 
is characterized by 

(14) U<t> {x) =mm{(f>(y) + p°(x -y) : y eV x }, x G O, 
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where V x = {y G dfl : ]y, x\ C Vt}. 

If (j) = 0, the Lax-Hopf function no is the distance from <90 with respect to 
the convex metric associated to p°, that is uq(x) = mm{p°(x — y) : y G 50}, 



if!! (see ]16j for a complete treatment of this subject). 
More generally, if (ft : d£l — > R satisfies 

(15) 0(x) - 0(y) < ^(x - y) Vx, y G 0J2, 
the Lax-Hopf function coincides with 

(16) u(x) := min{0(y) + p (x — y) : y G 90} , x G . 
Namely, from ([l^) it is easy to check that u = <fi = Us on dVL. Moreover, 



from the representation formula (14) it is clear that u < u^. On the other 
hand, let x G and let y G Oft be such that u(x) = 4>(y) + p°(x — y). Let 
z G \y, x\ C\dVL be such that }z,x{ C O. Then, by the representation formula 
(0) and by (|TJ) we have that 

u(x) < u^(a?) < <j)(z) + p°(x -z)< <f>(y) + p°(z-y) + p°(x - z) 

= 4>(y) + P°(x - y) = u(x), 

so that u(x) = u^(x). 

The following example shows that if we merely require (H4), then in 
general ^ u. 

Example 4.3. Given e G (0, tt) let us define the set 

O := {(pcos0,psin0) G R 2 : l<p<2, -tt < 9 < tt - e} , 

and let <f>: d£l — > R be defined by 4> = 8. If we choose K = Bi(0), then 
p°(£) = |£| for every £ G R 2 , and £(7) is the usual (Euclidean) length of 
a curve 7. A straightforward computation shows that <j> satisfies (H4). On 
the other hand, it is also clear that (j) does not satisfy (|l5"|). Namely, if we 
consider the two boundary points y\ = ( — 1, 0), y% = (cos(7r — e), sin(-7r — e)), 
then we have 

4>{V2) - <P{yi) = 2vr - e > 2sin(e/2) = \y 2 - y\\. 

In this case, the function u defined in (|lrj|) does not belong to since 
u(y2) < 0(yi) H - 1 2/1 — 2/2 1 < 0(2/2) and hence it does not achieve the boundary 
value at yi- 

We conclude this section with the following result due to S. Bianchini 
(see H and for its generalization to the case of nonsmooth constraint 
sets K). 

Theorem 4.4. There exists a weak solution Vf G L\(Q) to i/ie transport 
equation 

(17) -div(vfDp(Du l / > )) = f inn. 

As a corollary of Theorems [4.2| and [4.4| we conclude that (u^,Vf) is a 



solution to (10) 
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5. The minimum problem 

In this section we shall discuss the connection between the system of 
PDEs ( JTo| ) and the variational problem 



(18) min|— J fudx: u G X$ 

It is clear that, since / is non-negative and ua, is the maximal element in 



Xgj,, then Ugh is a solution to (18). We shall show that (|10| ) is the Euler- 
Lagrange necessary and sufficient condition for the minimum problem ( |l8|) . 
This property, which is the analogous to the duality in Mass Transport 
Theory, will be useful in the proofs of the uniqueness of the solutions of 

©• * 

As a preliminary step, we need to enrich the class of test functions allowed 
in ©. 

Lemma 5.1. Let u G C{£1) n W 1,0 °(Q), u = on dQ. Then there exists a 
sequence (tpj) C Cq°(Q) such that tpj — > u and Dipj — > Du a.e. in £1, and 
\\ipj\\oo < \\u\\oo, \\DipjWoo < 3||Z>u||oo for every j G N. 

Proof. We use a standard truncation argument. Let G G C 1 (R) such that 
< G(t) < \t\ for every t G M, G(t) = for |t| < 1, G(t) = \t\ for \t\ > 2, 
\G'(t)\ < 3 for every t G R. 

Let us define Uj := G(ju)/j, so that Uj — > u pointwise in £2. Moreover, 
Duj = G'(ju)Du, hence Duj — > Du a.e. in Q and \Duj\ < 3|Dn|. 

Let <p e be the standard family of mollifiers in R". Since Uj has compact 
support in Q, we can choose a sequence Ej \ 0, such that the sequence 
tpj := (p ej * Uj has the required properties. □ 

Corollary 5.2. C^°(Q) is dense in the space X := {u — w : u,w G X^} 

with respect to the weak* topology ofW 1,oc (Q). As a consequence, if(u,v) G 
X^ x L+(Q) satisfies fcH\), then 

(19) / v(Dp(Du), Dw - Du) = / f(w - u) \/w G X+. 

Jn Jn 

We are now in a position to prove that the Monge-Kantorovich system 
(fill) is the Euler-Lagrange condition for the minimum problem (18). 



Theorem 5.3. The minimum problem l\Tq ) and the system of PDEs (1C) 
are equivalent in the following sense. 

(i) u G Xfj, is a solution to fcH\) if and only if there exists v G Li_(Q) 
such that (u, v) is a solution to (j7~ 



(ii) Let u G Xa be a solution to fjlffi). Then (u,v) is a solution to (1C) 
if and only if (u<f,,v) is a solution to (]T~ 



L K (x) 



Proof. Let us denote by Lk the indicator function of the set K, that is 

Jo ifx€K, 
\ +oo if x G R n \ K , 
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so that 

F(u) := [ [Ik(Du) - fu]dx = - [ fudx, Vu£l f 
Jn Jn 

Since p is differentiable in M n \{0}, the subgradient of the indicator function 

can be explicitly computed, obtaining 

'{aDp{£): a > 0} if £ G dK, 
if £ K, 

{0} if feint If 

(see e.g. g Sect. 23]). 

Hence, if (u,v) is a solution to (|10|) , then v(x)D p(Du(x)) G 3Ik{Du(x)) 
for a.e. x G O, so that, for every u; G 

(21) F(w)-F(u)> / v(Dp(Du), Dw - Du) dx - / / (w - u) dx = , 



(20) 0/^(0 



where the last equality follows from Corollary |5.2j . This proves that u is a 
solution to fll8|) . 

Assume now that u G -X^ is a minimizer for i 7 , so that f(u — us) = 
a.e. in f2, due to the maximality of us in and the fact that / > 0. By 
Corollary 5.2 we can choose u — u^ as test function in the transport equation 
( |l7| ) solved by u^, getting 

= / v f {Dp(D U(j} ), Du- Du^dx = - / v f (l - {Dp(D U(j} ), Du)) dx . 
Jn Jn 

On the other hand, by Theorem ^]l](ii) and the fact that Du G K a.e. in S7, 

we have 

1 - (Dp(Du ( j > (x)), Du(x)) > 0, a.e. x G , 

1 - (Dp(D U(j) (x)), Du(x)) = Dp{D U4 ,{x)) = Dp(Du(x)), 

so that VfDp(Du ( p) = VfDp(Du) and uj(l — p(Du)) = a.e. in f2, that is 
(u,vj) is a solution of (|10D. This concludes the proof of (i). 

Let us prove (ii). The previous computation shows that if (u^,,v) is a 
solution of ([H]), and u G X$ is a solution of the minimum problem (18), 
then also (u,v) is a solution of (|T0|). Finally, let (u,u) be a solution to (10). 
Upon observing that vp(Du) = v a.e. in f2, and choosing u — u^ as test 
function in the weak formulation of the transport equation 

-dw(vDp(Du)) = f, in O 

we conclude that the opposite implication holds. □ 

6. Uniqueness of the ^-component 



Thanks to Theorem |5.3| (ii), the uniqueness of the v component of the so- 



lution of the system of PDEs 
Vf introduced in Theorem ^ 
equation 

(22) -dw(vDp(D U(j> )) = f. 



10] ) will follow from the fact that the function 
is the unique weak solution in L+(Cl) of the 
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We shall use an explicit representation of vt proved in || and, in order 
to explain it, we need to introduce some additional tools, related to the 
directions where the function u<f> has the maximal slope. 

Recalling the representation formula (|l4|) for u^, for every x S O let us 
define the projections of x by 

n(x) := {y G V x : U(f> (x) = <f>(y) + p°{x - y)} 

and, for every x G 0, let A(x) be the set of directions through x 

A(x) := { X ~ V : yeU(x)) , x G ft . 

Let I? C f2 be the set of points with multiple projections, that is 

(23) D := {x G SI : A (a;) is not a singleton}, 

and for every x G O \ D, let p(x) and d(x) denote the unique elements in 

II(x) and A(x) respectively. Finally, let r(x) be defined by 

(24) 

sup{i > 0; u^(x + sd(x))) = ii^(x) + s Vs G [0, t]}, x G O \ D , 
x £ D , 



r(x) : = 



and let J be the set 

(25) J := q(x), q(x) := x + 7~(x)d(x), 
where we understand that g(x) = x if x G -D. 

Definition 6.1. We shall call transport ray through x G any segment 

b,g(x)], pg n(x). 

It is clear from the definition that, if x G \ D, then there is a unique 
transport ray |p(x), q(x)J through x. On the other hand, if x G I?, any 
segment [p, x], with p G n(x), is a transport ray through x. 

The transport rays correspond to the segments where grows linearly 
with maximal slope, as stated in the following lemma. Hence in the sandpiles 
model, they correspond to the directions where the matter is allowed to run 
down, that is where the u-component could be nonzero. 

Lemma 6.2. For every if!! and y G n(x) we have 

(26) u+{z) = <P(y) + p\z -y) Vz G [y, xj . 
Moreover, if x G O \ D , then 

(27) u^z) = <j){p{x)) + p°(z - p(x)) \/z G [p(x), q(x)] . 

Proof. The identity (|26| ) follows from @, while (^?j) for x £ ft\ D follows 



from (|26|) when z G [p(x),x], and from the very definition of q(x) when 
ze{x,q(x)j. □ 

Corollary 6.3. Given x G Q\D, let {p(x), q(x)J be the transport ray through 
x. Then H(z) = {p(x)} for every z G ]p(x), g(x)[. 
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Proof. Prom (27) it is clear that p(x) G H(z) for every z G ]p(x), g(x)[. 
Assume by contradiction that there exists z G ]p(x),g(x)[ such that U(z) 
contains a point y ^ p(x). 

We claim that the three points y, z and q{x) cannot be aligned. Namely, 
we cannot have z G Jy, q(x){, for otherwise we would have y = p(x). On the 
other hand, we cannot have q(x) G indeed, in this case, by (^?j) we 

would have u^{q{x)) = 4>(y) + p°(q(x) — y), while, by (p6|), we know that 
u ( f ) (q{x)) = u ( p(z) + p°(q(x) — z). Putting together these information, we get 

u$(z) = 4>{y) + p°{z -y) = 4>(y) - p°(z - q{x)) + p°(q(x) - y) 
= u^z) + p°(q(x) - z)+ p°(z - q(x)) > u^z), 
a contradiction. 

Since the curve [y, z\ U \z, q(x)J, joining y to q(x) and lying in £1 except 
for its endpoints, is not a segment, by Lemma |2.4| it cannot be a geodesic, 
so that 

U4>{q{x)) < (f>(y) + p°(z -y)+ p°(q(x) - z). 
Taking into account ( |27| ) we get 

u^z) = u^qix)) - p°{q(x) - z) < 4>(y) + p°{z - y), 
a contradiction with y G n(z). □ 

At this stage we are dealing with three meaningful sets of singular points 
associated to u^: 

S, the sets of points in O where is not differentiable; 
D, the sets of points in with multiple projections; 
J, the set of endpoints in £1 of the transport rays. 

The set J plays a central role in optimal transport theory. As a matter 
of fact, the disintegration of the Lebesgue measure in the reference set f2 
along transport rays, can be achieved only if J has Lebesgue measure zero, 
and this is one of the main tools of the theory. 

If and K are sufficiently smooth (e.g. both of class C 2 , and K with 
strictly positive curvatures), one can prove that D coincides with E (see, e.g. 



[14, 16]), so that, by Rademacher's Theorem, D has zero Lebesgue measure. 
Moreover, using the methods of Differential Geometry, it is possible to prove 
that J = £ and to characterize J \ S in terms of optimal focal points. As 
a consequence one can prove that J has zero Lebesgue measure and give 



detailed rectifiability results (see, e.g., [16, [21|). 

If Q, and K are not regular, the set X may have positive Lebesgue measure 
(see [^] for an example where G C 1,Qf , a < 1, and K is the unit ball), 
and J can be a proper subset of E. 

We recall here some known properties of the singular sets, adding some 
remarks in order to complete the description of the relationships between 
these sets. The known results, proved in ||, are collected in the following 
proposition. 

Proposition 6.4. The sets E, D and J satisfy the following properties. 
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Figure 2. The sets ft and K v of Example 6.5 



(i) D CE, More precisely, if is differentiable at x £ ft, then x has 
a unique projection and A(x) = {Dp(Du < j } (x))} . 

(ii) The set J has zero Lebesgue measure. 

(iii) D is a (n — 1) rectifiable set. 

Remark 6.5. We recall that D = £ if K is strictly convex. Proposition 
6.4(i) states that D C S, in general. The following example shows that it 
can happen that D ^ S. Let us consider the set 

2 



(x,y) EM 2 : Mx| + ^ | + //" 1 1 



and let K be its polar set. The set K° is strictly convex but it is not of 
class C , hence the set X is of class C 1 but it is not strictly convex (see 
Theorem 2.1). Let O be the hexagon 

n = {{x,y) £R 2 : \y\ <2-\x\, \y\ < l} 

(see Figure §). The Lax-Hopf function uq associated to K and to the 
boundary datum <f> = is the distance function to the boundary of f2 with 
respect to the convex metric p°, which can be explicitly computed as 

u (x) = /V2(l-|y|), if (x,y) £nn{\x\ < 1}, 

U ° W \V2(2-\x\-\y\), i£(x,y)eSln{\x\>l}, 

so that S \ D = {(x,y) £ ft: |x| = 1, y ^ 0}. 



Remark 6.6. Since D is a subset of J, then by Proposition |6.4j (ii) it has 
Lebesgue measure zero. Notice that D can be a proper subset of J even in 
the regular case. For example, let ft be the ellipsis {(x,y) £ M 2 ; x 2 /a? + 
y 2 /b 2 < 1}, where < b < a, and let K = Bi(0). The points P = 
((a 2 — 6 2 )/a,0) and Q = (—(a 2 — b 2 )/a,0) are the centers of curvature of 
9ft at (a, 0) and (—a, 0) respectively. It can be checked that J = {Q, PJ 
whereas D = }Q,P{. 

Moreover, since J need not be closed (see ||, Remark 5.2]), J 7^ D in 
general. The following proposition shows that JC D. 
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Proposition 6.7. For every x G Q \ D we have that t(x) > 0, where r is 



the function defined in (24). 



Proof. Let us fix x G VL \ D. We can assume without loss of generality that 
x = 0, and d(0) = —Xe n , A > 0. Since the map z h-» d(z) is continuous 
in the open set O \ D (see |], Prop. 3.2]), there exists r > such that 
B r (0) cn\D, and (d(z), e„) < -A/2 for every z ejB r (0). 

Given 5 G (0,r/2), let us define the function jB r (0) -> R n by 

(28) ^) == z + f 2 f~^' , 



(d(z), e n ) 

so that (ift(z), e n ) = 25. Hence ifi maps B r (0) into the hyperplane H := 
{x G R n : (x, e n ) = 25}. The continuity of d(z) in B r (0) ensures that we 
can choose 5 in such a way that ip(z) G B r (0) for every z G Bs(0). 

Summarizing, ^: B r (0) — > H is a continuous map, ift(Bs(0)) C B r (0), 
and ■0( 2; ) belongs to the ray {z — td(z) : t > 0}. Moreover, for every 
2 € -B<5(0), we have that {ip(z), z\ C B r (0) C fi, hence ij){z) G and, 



by Corollary 6.3, we infer that d(z) = d{ip(z)). 

Since iff maps dBs(0) into the (n — l)-dimensional hyperplane H, by the 
Borsuk-Ulam lemma (see [ pO"| , Cor. 4.2]) there exists w G dB$(0) such that 
t/>(u;) = ip(—w). Furthermore <i(w) = d(tp(w)) = d(ip(—w)) = d(—w), so 
that the equality ip(w) = ip(—w) implies that 

(w, e n ) 
{d(w), e n ) 

i.e., w is parallel to d(w). But the transport ray through w and — w contains 
the origin, hence d(w) = d(—w) = d(0) and r(0) > 5. □ 

For the reader's convenience we collect here the results concerning the 
relationships between the singular sets. 

Corollary 6.8. The sets E, D and J have zero Lebesgue measure. More- 
over, D C S and D C J C D, with possibly strict inclusions. 

We are now in a position to briefly recall the explicit representation of 
the function v j appearing in ([17]) , one of the main results in . 

Theorem 6.9. There exists a positive function a : f2 — > (0, +oo), absolutely 
continuous along almost every transport ray, such that an explicit solution 
Vf G L+(Q) of ( fj^ ) is given by 

(29 Vf (x)= f(x + td(x))— , N V " dt, a.e. xen. 

Jo ' 

More precisely, for a.e. x G Q \ D the function Vf is locally absolutely 
continuous along the ray 1 1— > x + td(x), t G [0,r(x)) ; and satisfies 



(30) 



— [vf(x + tdix)) a(x + id(x))] = — /(x + tdix)) aix + tdix)), 
dt 

lim u/(x + a(x + id(a;)) = 0. 

t-tr(x)~ 
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Proof. See ||, Theorems 5.7 and 7.1. □ 
Remark 6.10. The reader may want to compare the representation for- 



mula (|29|) with the one proved in the regular homogeneous case (see [1C, 



11, p.q|), where a(x) can be explicitly written in terms of the distance from 
the boundary uo(x) and of the anisotropic curvatures m, . . . , ft n -i of dCl at 
p(x). More precisely, one has 

n-1 

a ( x ) = II — u o( x ) K j(p( x ))] > X £ fl\D . 
3=1 

We shall find a necessary and sufficient condition for having that Vf is the 
unique solution to (|22|). The following example shows why the uniqueness 
of the ^-component may fail. 

Example 6.11. Let O = (0,1) x (0,1) C M 2 , K = Bi(0), and let 0(x,y) = y 
on dQ, so that iu(x, y) = y in Then every function of the form v(x, y) = 
l — y+c(x), c(x) > for x G [0, 1], is a non-negative solution of the equation 
— div(v Du^) = 1. Notice that in this case Vf(x, y) = 1— y in U, so that (p2|) 
has infinitely many solutions of the form v = Vf + v, with div(v Du^) = 0. 
The existence of non-trivial variations of vj is due to the fact that Q is 
covered by transport rays with both endpoints on dQ. 

Now it should be clear that, in order to discuss the uniqueness of the 
t;-component of the solutions to (]22|), we need to introduce the set 



where the set E is defined by 

E := {{p, q)Gdnx On; p^q, jp, q{ C cp(q) = <f>(p) + p°(q - p)}. 
In Example |6.11 , E = {((x,0), (x, 1)), x G (0, 1)} and T = Jl. In general 



T is the union of the transport rays ending on dQ.. On this set every function 
in is forced to have maximal slope and hence to coincide with ua. 

Lemma 6.12. Let (p,q) G E be given. Then the following hold. 

i) II(x) = {p} and q(x) = q for every x £ Jp, q\. 

ii) If u G X^, then u{x) = <p(p) + p°(x — p) for every x G [p, g]. As a 
consequence every u G coincides with on the set T. 

Proof. Let us start by proving that p G II(x) for every x G Jp, q{. Namely 

U^ix) < <f>{p) + p°(x -p) = <j){q) - p°(q -p) + p°(x - p) 

- p°(q -x)< u^x) , 



where the last inequality follows from Lemma ^2. Hence uAx) = 4>(p) + 
p°(x — p), so that p G II(x). 

Actually we can prove that II(x) = {p} for every x G Jp, q{ using the same 
argument of the proof of Corollary |6.3|: if there exists p\ ^ p, p\ G II(x), 
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then the three points pi, x and q cannot be aligned, so that 
4>(q) - 4>(pi) < p°(q -x)+ p°(x - p^ 

and 

M x ) = ^(pi) + p°( x -Pi) > 4>(q) - p°(q -x) = 4>(p) + p°(x -p) = u^x), 

a contradiction. 

Since q G <9f2 and, for every x G ]p, </[, H(x) = {p}, it is clear that 
q(x) = q. 



we 



Finally, let u G and x G [p, (/] be given. Then, using Lemma 2^2 
get 

(f>(p) + -p) = 4>{q) - p°(q -x)< u(x) < (f>(p) + p°(x - p) , 
and (ii) follows. □ 



Theorem 6.13. The function vj defined in $j9Q is the unique solution of 
Mm if and only if the set T has zero Lebesgue measure. 



Proof. The uniqueness of v f when T has zero Lebesgue measure is proved in 
||, Prop. 7.3]. The proof is based on the fact that every weak solution v G 
L l + {Q) to ( |22| ) is locally absolutely continuous along almost every transport 
ray, and satisfies 

— [v(x + td{x)) a(x + td(x))\ = -f(x + td(x)) a(x + td{x)) . 

Moreover, if the endpoint q{x) of the transport ray through x belongs to ft, 
the v ■ a satisfies the initial condition 

lim v j(x + td(x)) a(x + td(x)) = 



so that, by Theorem |6.9| , v = v t along that ray. 

Let us now assume that T has positive Lebesgue measure. It is clear that 
the function 

w^(x) := mm{-(p(y) + p°(y - x) : y G V x }, xGO, 

is the Lax-Hopf function corresponding to the geometry induced by the 
convex set —K and the boundary datum —<j). Assume that and are 
differentiable at a point x G }p, q{, with (p,q) G E. From Proposition 6A 
we have that Du^(x) = (q — p)j p°(q — p). The same argument, used with 
—<p as boundary datum and the geometry induced by —K, shows that 

Dp-(Dw^(x)) = p0 ^ p ^ g) = - p o\ q I p) = Dp(Du^(x)), 

where p- and p®_ denote respectively the gauge functions of the convex sets 
-K and -K°. 

Let us consider the indicator function of the set T, 



g(x) :-- 



1, x G T, 

o, x g n \ T. 
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By Theorem |4.4| , there exist two nonnegative functions Vg , v g such that 
- div (v g ^Dp(Du ( j > )) = g , - div(v g Dp- (Dw^)) = g, in Q . 



Looking at the explicit representations of and v as in (p9|), we easily 



get 

v g h =v~ = a.e. in \ T , v J > 0, v~ > a.e. in T , 
so that, for every p G CQ°(f2), we obtain 

(«a + Vg)(Dp(D U4> ), Dp) dx = J^(v+ + v g )(Dp(Du^), Dp) dx 

v~ g h {Dp{Du 4> ), Dp)dx - / v~ (Dp-(Dw (j) ), Dp) dx 
t Jt 



gpdx — I gpdx = . 

Hence for every A > the function v\ = Vf + + is a non-negative 
weak solution to (22). □ 

Remark 6.14. It can be easily checked that if the datum <p satisfies (H4) 
with strict inequality holding for every x ^ y, then the Lebesgue measure 
of T is zero, so that Vf is the unique w-component allowed in ([R 



Remark 6.15. If the Lebesgue measure of Q\T is zero, then by Lemma 5.12 
we have that Xs = {u^}. On the other hand, it can be easily checked that 
the converse implication is also true. In ||] (see also H) it is proved that 
either / {u^} or there exists a weak solution v G Lf oc (Q) to 

(31) div(v D U(t> ) = 

such that v > a.e. in f2. In the proof of Theorem 3,13| we have shown 



that, as soon as T has positive Lebesgue measure, it is possible to construct 
a solution v G L+(£l) to (|3ll ) such that v > a.e. on T. 

7. Uniqueness of the u-component 



In Theorem |5.3| we have shown that the set of admissible -u-components 
of the solutions to (|l0|) coincides with the set of solutions to the minimum 



problem (18) 



In this section we shall construct, using Lemma 3.1, the minimal solution 
to (|T8|), and we shall give a necessary and sufficient condition in order to 
have that this function equals the maximal solution in f2. 

Let us define the function uf. fi-^Rby 

(32) u f (x) : = sup{u^(z) - £(7) : z G dil U supp(/), 7 G L^}, 

where supp(/) is the complement in Q of the union of all relatively open 
subsets A C Q such that / = a.e. in A. 

Proposition 7.1. The function Uf is characterized by 

(33) u f( x ) = max{tt^(z) — p°(z — x) : z G W^}, 

where W x := {y G <90 U supp(/) : ]y, x[ C O \ supp(/)}, and satisfies 
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(i) u f E X^; 

(ii) Uf[x) = u ( p(x) for every x E supp(/) U T; 

(iii) given u E X^, we have u = on supp(/) if and only ifuf < u < 
in O. 

Proof. The function Uf is the function defined in (|7|) corresponding to 
D = Q, C = dft U supp(/), and ip = on C. Clearly the function ip 
satisfies (@) on C, since u$ satisfies that condition on Q. Then the stated 



properties follow from Lemmas 3.1, 3.3 and 6.12. □ 



Theorem 7.2. A function u E X^ is a solution to (18) if and only if 



< u < u^. Moreover, the function uj coincides with in Q (and hence 



is the unique solution to (18)) if and only if J HQ C supp(/). 



Proof. The first assertion follows from Proposition |7.1| (iii) and from the 
fact that u E X^ is a solution to (pjj|) if and only if u = on supp(/). 
Let us prove the uniqueness result. Assume that Jflfl C supp(/), so that 



J E 9$7Usupp(/) and, from Proposition |7. l|(ii) , Uf = on J. Let x E fi\ J 



be given, and let q(x) E J be the endpoint of the ray through x. We have 

u^x) = u<j,(q(x))- p°(q(x)-x) = u f (q(x)) - p° (q(x) - x) < Uf(x) < u$(x) , 

and hence uj(x) = u^x). 

Assume now by contradiction that Uf = in Q and that there exists 
xq E J n ft such that xo supp(/). Let z E W xo and y E V xo be such that 

(34) Uf(x ) = u^z) - p°(z - x ), 

(35) u^xq) = u<t,{y) + p°(x - y). 

Notice that y / z, otherwise we get —p°(y — xo) = p°(xo — y). Hence xq, y 
and z are three distinct points. Moreover we have that xq E Jy, z\, otherwise 



by Lemma 2.4 we should have 

u^z) - u^y) < p°(z - x ) + p°(x - y) 

and 

u^xq) = u^z) - p°(z - x ) < u^y) + p°(x - y) = u^xo), 

a contradiction. 

Finally, the fact that xo E ]y, z\ implies that 

u^z) = u^xq) + p°(z - x ) = u^(y) + p°(x - y) + p°(z - x ) 

= u ct> (y) + p (z-y), 

hence ]y, z\ is a transport ray through xq, so that xo 2" J ■ □ 
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